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Abstract — This paper presents two methods for approximating 
the performance of coded multicarrier systems operating over 
frequency-selective, quasi-static fading channels with non-ideal 
interleaving. The first method is based on approximating the 
performance of the system over each realization of the channel, 
and is suitable for obtaining the outage performance of this type 
of system. The second method is based on knowledge of the 
correlation matrix of the frequency-domain channel gains and 
can be used to directly obtain the average performance. Both of 
the methods are applicable for convolutionally-coded interleaved 
systems employing Quadrature Amplitude Modulation (QAM). 
As examples, both methods are used to study the performance 
of the Multiband Orthogonal Frequency Division Multiplexing 
(OFDM) proposal for high data-rate Ultra-Wideband (UWB) 
communication. 

I. Introduction 

Multicarrier communications systems based on Orthogonal 
Frequency Division Multiplexing (OFDM) have gained inter- 
est from the communications community in recent years, as 
evidenced by wireless standards such as 802.11 a/g and high- 
rate Ultra-Wideband (UWB) [1], [2]. In many situations, the 
wireless channel for these systems can be assumed very slowly 
time-varying relative to the transmission rate of the device, 
and can be approximated as quasi-static for the duration 
of one or more packet transmissions. As well, since most 
OFDM devices use a relatively large bandwidth, the channel is 
frequency-selective. This motivates an interest in the analysis 
of the performance of coded multicarrier systems when trans- 
mitting over quasi-static frequency-selective fading channels. 
Throughout this paper, we will assume that the OFDM system 
is designed such that the cyclic prefix is always longer than the 
channel impulse response. Thus, we can consider the channel 
in the frequency domain, with one channel gain per subcarrier 
and a correlation matrix representing the correlation between 
the gains of different subcarriers. 

When considering the performance of systems over quasi- 
static channels, the concept of outage is often used. That is, 
since (in a packet-based transmission system) each packet will 
be transmitted over only one realization of the quasi-static 
channel, we consider some channel realizations to be in outage 
(that is, they do not support the required data rate). The worst- 
case performance of the non-outage cases is then studied, 
providing information about what minimum performance can 
be expected of the system given a certain allowable outage 
rate. 



It should be emphasized that classical bit error rate (BER) 
analysis techniques for coded systems [3], [4] are not appli- 
cable in this setting because the channel is (a) non-ideally 
interleaved (resulting in non-zero correlation between adjacent 
coded bits), and (b) quasi-static (which limits the number of 
distinct channel gains to the number of OFDM tones). 

Related Work: In [5], Malkamaki and Leib consider the per- 
formance of convolutional codes with non-ideal interleaving 
over block fading channels. They make use of the generalized 
transfer function (GTF) [6] in order to obtain the pairwise 
error probability in the case when the channel gains for each 
transmitted symbol are not unique (due to a small number 
of fading gains as a result of the block fading model). 
This approach is similar in some ways to that presented in 
Section UlI-CI The major difference is that our method does 
not require the GTF of the code, which may be difficult to 
obtain [7]. Instead, we apply the concept of error vectors, 
introduced in Section UlI-AI 

In [8], Veeravalli presents a general approach for analyzing 
the performance of uncoded diversity-combining systems on 
correlated fading channels. In particular, the probability of 
error is expressed (for Rayleigh-fading channels) as a function 
of the complex covariance matrix of the channel gains. In 
Section UlI-DI we extend this approach to handle coded mul- 
ticarrier systems employing arbitrary Quadrature Amplitude 
Modulation (QAM) modulation schemes. 

Organization: Section [II] introduces the system model con- 
sidered in this paper. In Section|lIl] we present two approaches 
to performance analysis: a per-realization based approach 
for system performance evaluation which readily leads to 
estimation of outage-based performance measures, and an 
approach, based on the channel correlation matrix, which can 
be used to obtain the average system performance. Numerical 
results illustrating the effectiveness of both proposed methods 
are presented in Section [TV] Finally, Section IVl concludes the 
paper 

II. System Model 

Throughout this paper we consider an A^-tone OFDM 
system, using A/-ary QAM (M-QAM) on each subcarrier The 
system employs a punctured convolutional code of rate Re- 

We assume that the transmitter selects a vector of 
RcN\og2M random message bits for transmission, denoted 
by b = [bi b2 ■ ■ ■bR^Niog^Ai]'^ (where [-J^ denotes vector 



transposition). The vectors c and of length ~ N\og2 M 
represent the bits after encoding/puncturing and after interleav- 
ing, respectively. The bits are then modulated using M- 
QAM on each subcarrier. The resulting TV modulated symbols 
are denoted by the vector x ^[xi X2 ... . 

The symbols x are transmitted through a quasi-static 
fading channel with frequency-domain channel gains h = 
[hi /i2 ... hisi]- Writing H = diag(ft,), where diag(ft.) 
denotes a matrix with the elements of h on the main diagonal, 
we can express the received symbols r as 

r^^/E'sHx + n, (1) 

where n is a vector of independent complex additive white 
Gaussian noise (AWGN) variables with variance J\fo and Es is 
the energy per modulated symbol. The energy per information 
bit is Eb^ EJiR^log^M). 

We assume perfect timing and frequency synchronization. 
The receiver employs a soft-output detector followed by 
deinterleaving, depuncturing, and Viterbi decoding, resulting 
in an estimate b = [foi 62 • • ■ &_Rciviog2 a/]"^ of the original 
transmitted information bits. 

III. Two Performance Analysis Methods 

In this section, we present two methods for approximating 
the performance of coded multicarrier systems operating over 
frequency-selective, quasi-static fading channels. The first 
method is based on approximating the performance of the 
system over individual channel realizations, and is suitable 
for obtaining the outage performance. The second method is 
based on knowledge of the correlation matrix of the frequency- 
domain channel gains and can be used to directly obtain the 
average performance. Both methods are based on considering 
the set of error vectors, introduced below. 

One major problem in the analysis of A/-QAM modulation 
schemes with A/ > 4 is that the probability of error for a 
given bit depends on the whole transmitted symbol (i.e. it also 
depends on the other bits in the symbol). For this reason, for 
the combination of convolutional coding and A/-QAM it is not 
sufficient to adopt the classical approach of considering devia- 
tions from the all-zero codeword. In theory, one must average 
over all possible choices for c. Since this is computationally 
intractable, we simply assume the transmitted information bits 
b (and hence x) are chosen randomly. For 2-QAM and 4- 
QAM (where the joint hnearity of code and modulator is 
maintained) this is exactly equivalent to considering an all- 
zero codeword. In the case of M > 4, we have verified for 
the two analysis methods proposed below that the results are 
practically invariant to the choice of b. 

A. Set of Error Vectors 

Let £ be the set of all L vectors Cj (1<7<L) of code 
output (after puncturing) associated with input sequences with 
Hamming weight less than Wmax, i-e., £ = {ei, 62, . . . , e/,}. 
Let Ij and aj be the length of ej and the number of information 
bit errors associated with ej, respectively. As with standard 



union-bound techniques for convolutional codes [3], the low- 
weight terms will dominate the error probability. Hence it 
is sufficient to choose a small Wmax — for example, the 
punctured code of rate 1/2 considered in Section Hvl has a 
free distance of 9, and choosing i(;,nax = 14 (resulting in a set 
of i = 242 error vectors of maximum length I ~ 60) provides 
results virtually identical to those obtained using larger Wmax 
values. 

In the following, we refer to ej as an "error vector" 
and to £ as the set of error vectors. Note that £ can be 
straightforwardly obtained from the transfer function of the 
code, without resorting to the GTF approach as in [5]. It is 
also independent of the number of distinct channel gains (or 
the number of blocks in the context of [5]). 

B. Pairwise Error Probability 

We consider error events starting in a given position i of 
the codeword {1 < i < Lc). We consider each error vector ej 
for 1 < j < L, and form the full error codeword 

q,^j = [0 ... ... 0]^ (2) 

of length Lc by padding Bj with zeros on both sides as 
indicated above. Given the error codeword ^ and given that 
codeword c is transmitted, the competing codeword is given 
by 

V, J = c e q,_j- (3) 

where © denotes XOR. Letting Zij be the vector of QAM 
symbols associated with (the interleaved version of fi.j), 
the pairwise error probability (PEP) for the jth error vector 
starting in the ith position is 

PEP,,j^Qy^J\H{x-z,M^^ . (4) 

where Q(-) is the Gaussian Q function [3]. 

C. Per- Realization Performance Analysis ("Method I") 

In this section, we obtain an approximation of the BER for a 
particular channel realization H = diag(/i), which we denote 
as P{H). This method leads naturally to the analysis of the 
outage performance, which is a useful measure for schemes 
operating over quasi-static channels. For simplicity, we shall 
refer to this as Method I in the remainder of this paper. 

The PEP for an error vector (1 < j < E) with the error 
event starting in a position « (1 < i < Lc} is given by Q. 
The corresponding bit error rate for this event is given by 

= •PEP,j(/J) . (5) 

Summing over all L error vectors, we obtain an approximation 
of the BER for the z"^ starting position as 

L 

P,(ff) = ^a, •PEP,j(ff) . (6) 

We note that (|6|l can be seen as a standard truncated union 
bound for convolutional codes (i.e., it is a sum over all error 



TABLE I 

Pseudocode for Method I. Final BER is P (for given H). 
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TABLE II 

Pseudocode for Method II. Final BER is P. 



P := 

for i :^ 1 to Lc do 
for j := 1 to L 
form q^j as per 
form Vij as per 

form ai'i j 





frt 

h.'ij and Rh'h' 



Compute D := diag(a;' — z'i.j) 
Compute Rgg := DR^'h'D^ 

Compute to obtain P; 

P := P + P,,, 
= P / 



events of Hamming weight less than ojmax)- We also note that 
we can tighten this bound by limiting Pi{H) to a maximum 
value of 1/2 before averaging over starting positions [5]. 
Finally, since all starting positions are equally likely, the BER 
P{H) can be written as 



P{H) 



min 

1=1 



1 ^ 



(7) 



Table U contains pseudocode to calculate P{H) according 
to 0. This method readily leads to the consideration of 
the outage performance. We evaluate (0 for many channel 
realizations. For a given X% outage rate, the worst-performing 
(100-X)% of realizations are considered in outage (i.e., they 
are not capable of supporting the required data rate), and the 
worst-case performance of the non-outage cases is shown. This 
provides information about the minimum performance that can 
be expected of the system given the X% outage rate. 

We should note here that does not take into account 
edge effects at the end of a coded OFDM symbol (i.e., near 
the end of a symbol not all error vectors represent valid errors 
because the codeword may end before the error vector does). 
It is possible to modify Q in order to disregard error events 
that do not "fit" within the symbol length. However, such 
modification makes a minimal difference in the final BER of 
the overall system. 

D. Averaging Approach ("Method U") 

In this section, we propose a method, based on knowledge 
of the frequency-domain channel correlation matrix, which 
can be used directly in order to obtain the average BER 
performance of coded multicarrier systems. For simplicity, we 
refer to this as Method II in the remainder of this paper. 

For this method we will explicitly assume that the elements 
of h are Rayleigh-distributed and have known correlation 
matrix T,hh (in practice, H^hh can be obtained from actual 
channel measurements, or can be numerically estimated by 
measuring the correlation in many channel reaUzations from 
a given model). 

Note that for an error vector starting in position i, only 
the non-zero terms of {x — Zij) contribute to the PER Let x' , 
z'ij and H' = diag(/i,') represent the transmitted symbols, 
received symbols, and channel gains corresponding to the non- 
zero distances of (x — Zij), respectively, and form Rh'h' by 
extracting the elements from S^h which correspond to h'. 



Let D = diag(a;' 
distances. Now g 



H'{x' - 

E{g) = 
E(gg") 



be the diagonal matrix of non-zero 
z', ,) = Dh'. We have 



Ra 



DRh'h'D 



H 



where £'(•) denotes expectation and (•)^ denotes Hermitian 
transpose. The distribution of g is zero-mean complex Gaus- 
sian with covariance matrix Rgg- Following [8, Eq. (7)], we 
can write the bit error probability for the j*'' error vector 
starting in the i'^ position as 



Pr 



tt/2 - 



det 



/ EgR, 



gg 



de 



(8) 



\Afo sin^ e 

Note that (|8} can be written in closed form in terms of the 
eigenvalues of Rgg, as in [8]. Alternatively, (|8j can simply be 
solved directly via standard numerical integration techniques. 
The latter approach is favorable in the case we are considering, 
since Rh'h' (and thus Rgg) is in general different for each 
necessitating a new eigendecomposition. 
Now, summing over all L error vectors, the BER for the i*'^ 
starting position can be written as 



P 



L 



(9) 



Finally, since all starting positions are equally likely to be 
used, the average BER P can be written as 



Le La L 

1=1 j = l 



(10) 



Table im contains pseudocode to calculate P according to JlOt . 

IV. Numerical Results 

In this section, we present numerical results illustrating the 
performance analysis methods presented in Section HiH 

A. System and Channel Model 

As an example OFDM system, we have chosen the Multi- 
band OFDM system (MB-OFDM) proposed for IEEE 802.15 
TG3a High Data Rate UWB [2]. The MB-OFDM system uses 
128 tones and operates by hopping over 3 sub-bands (one 
hop per OFDM symbol) in a predetermined pattern. We will 
assume that hopping pattern 1 of [2] is used (i.e. the sub-bands 



are hopped in order). As a result we can consider MB -OFDM 
as an equivalent 384 tone OFDM system. After disregarding 
pilot, guard, and other reserved subcarriers, we have N ~ 300 
data-carrying tones. 

Channel coding in the proposed standard consists of a 
punctured maximum free distance rate 1/3 constraint length 
7 convolutional encoder and a multi-stage block-based in- 
terleaver (see [2] for details). After modulation, modulated 
symbols are optionally repeated in time (in two consecutive 
OFDM symbols) and/or frequency (two tones within the 
same OFDM symbol), reducing the effective code rate by a 
factor of 2 or 4 and providing an additional spreading gain 
for low data rate modes. In the framework of our analysis, 
we can equivalently consider this time/frequency spreading 
as a lower-rate convolutional code with repeated generator 
polynomials. In the proposed standard, the interleaved coded 
bits are mapped to quaternary phase-shift keying (QPSK) 
symbols using Gray labeling. As an extension for higher data 
rates, we also consider Gray-labeled 16-QAM. 

For a meaningful performance analysis of the MB-OFDM 
proposal, we consider the channel model developed under 
IEEE 802.15 for UWB systems [9]. The channel impulse re- 
sponse is a based on a modified Saleh-Valenzuela model [10]. 
As well, the entire impulse response undergoes an "outer" log- 
normal shadowing. The channel impulse response is assumed 
time invariant during the transmission period of (at least) one 
packet (see [9] for a detailed description). We consider the 
UWB channel parameter set referred to as CMl [9]. 

The OFDM transmit signal experiences a frequency non- 
selective fading channel with fading along the frequency 
axis. Thereby, the outer lognormal shadowing term mentioned 
above is irrelevant for the fading characteristics as it affects 
all tones equally. Denoting the lognormal term by G, we 
obtain the corresponding normalized frequency-domain fading 
coefficients as 

hr = h/G. (11) 

The elements of /i" are well-approximated as zero-mean 
complex Gaussian random variables [11]. This allows us to 
apply the analysis of Section IllI-DI to the UWB channel 
without lognormal shadowing, and then average over the 
lognormal shadowing distribution in order to obtain the final 
system performance over the composite UWB channel. We 
note that this is only relevant for Method II — for Method I 
the distribution of h is not important. 

B. Realization-based Approach 

In FigureQlwe present the 10% outage BER as a function of 
10 logj^Q {Eb /Ao) (the signal-to-noise ratio per information bit) 
obtained using Method I (lines), as well as simulation results 
(markers) for different code rates and modulation schemes 
using a set of 100 UWB CMl channel realizations with 
lognormal shadowing. The 10% outage BER is a common 
performance measure in UWB systems, cf. e.g. [2]. We can see 
that Method I is able to accurately predict the BER for QPSK 
and 16-QAM modulation schemes with a variety of different 
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Fig. 1. 10% Outage BER vs. 10 logm (Ej/Ao) from Method I (lines) and 
simulation results (markers) for different code rates and modulation schemes. 
UWB CMl channel with lognormal shadowing. Code rates 1/4 and 1/8 
include repetition. 

code rates, with a maximum error of less than 0.5 dB. It is also 
important to note that obtaining the Method I result requires 
significantly less computation than is required to obtain the 
simulation results for all 100 UWB channel realizations. We 
should also note that it is not clear that even 100 channel 
realizations (standard for UWB system performance analysis) 
is sufficient to accurately capture the true system performance 
(see Section HV-DI for further discussion of this issue). 

C. Correlation-based Approach 

In Figure |2l we present the average BER as a function 
of 101ogiQ(£'6/7Vo) for different codes rates with 16-QAM 
modulation over UWB channel CMl. The solid lines are 
the average BER obtained directly from the method of Sec- 
tion Illl-DI while the dashed lines are comparisons obtained 
from applying the method of Section UlI-CI to 1000 channel 
realizations and averaging the results. We show curves both 
with and without the "outer" lognormal shadowing of the 
UWB channel model. As can be seen from Figure |2j when 
the lognormal shadowing is neglected, the results obtained 
from the two methods show a close agreement. The difference 
between the two curves with lognormal shadowing is due to 
a slackening of the average bound which occurs when the 
results of Method II are averaged numerically over the "outer" 
lognormal distribution. 

D. A Note of Caution for System Designers 

Due to the high complexity of simulating system per- 
formance over quasi-static channels (which necessitates a 
large number of simulations over many channel realizations), 
system designers are tempted to use a relatively small set 
of (for example) 100 channel realizations to estimate system 
performance. Figure |3j a) (dashed lines) indicates the average 
BER with respect to 101ogjQ(i?b/Ao) for 4 different sets of 
100 UWB CMl channel realizations, obtained via Method I. 
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Fig. 2. Average BER vs. 10 logjQ(£'i,/A/'o). Solid lines: correlation-based 
approach (Method II). Dashed lines: realization-based approach (Method I) 
for comparison. UWB CMl channel, with and without lognormal shadowing 
("LN" and "no LN", respectively). 16-QAM, code rates 1/2 and 3/4. 

For comparison, the average performance obtained via Method 
II is also shown (solid line). We can see that the average 
system performance obtained using sets of only 100 channel 
realizations depends greatly on the specific realizations which 
are chosen. Similarly, Figure |3jb) illustrates the 10% outage 
BER with respect to 10 logjo(£'b/A/o) for 4 different sets of 
100 UWB CMl channel realizations, obtained via Method I. 
For comparison the 10% outage BER obtained using a set of 
1000 realizations is also shown. We see that the outage BER 
curves, while less variable than the average BER curves in 
Figure |3la), are still quite dependent on the selected channel 
realization set. 

Based on the results above, it seems that performance 
evaluation for systems operating on quasi-static channels using 
only small numbers of channel realizations may be prone to 
inaccurate results. This is one of the main strengths of Method 
I presented in Section IIII-CI the performance can easily be 
evaluated over any (large) number of channel realizations 
without resorting to lengthly simulations. 

V. Conclusions 

In this paper, we have presented two methods for eval- 
uating the performance of convolutionally-coded multicar- 
rier systems employing QAM and operating over frequency- 
selective, quasi-static, non-ideally interleaved fading chan- 
nels. The realization-based method ("Method I") presented 
in Section IIII-CI estimates the system performance over each 
realization of a channel with an arbitrary fading distribution, 
and is suitable for evaluating the outage performance of 
systems. The method presented in Section lTlI-DI CMethod 11"), 
based on knowledge of the correlation matrix of the frequency- 
domain channel gains, allows for direct calculation of the 
average system performance over the ensemble of quasi-static 
Rayleigh fading channel realizations. Both methods provide 
accurate estimates of the system performance as demonstrated 




(a) 101ogi(,(B6/A/'o) . (b) 101ogio(Ei,/AAo) . 

Fig. 3. (a) Average BER vs. 10 logi,)(Ei,/Ao) for four different sets of 100 
channels using Method I (dashed lines). For comparison: average performance 
from Method II (solid line), (b) 10% Outage BER vs. 10 logioCSfc/A^o) 
for four different sets of 100 channels using Method I (dashed lines). For 
comparison: 10% Outage BER for a set of 1000 channels (solid line). Code 
rate 1/2, 16-QAM, UWB CMl channel with lognormal shadowing. 

by the numerical results of Section IIVI and are both much 
simpler to evaluate than attempting system simulations over 
large sets of channel realizations. 
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